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^ ■ Abstract 

ff) ■ Cells sense the geometry and stiffness of their adhesive environment by active contractility. 

For strong adhesion to flat substrates, two-dimensional contractile network models can be used 
to understand how force is distributed throughout the cell. Here we compare the shape and 
force distribution for different variants of such network models. In contrast to Hookean networks, 
,-0 , cable networks reflect the asymmetric response of biopolymers to tension versus compression. For 

passive networks, contractility is modeled by a reduced resting length of the mechanical links. 
In actively contracting networks, a constant force couple is introduced into each link in order to 
model contraction by molecular motors. If combined with fixed adhesion sites, all network models 
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lead to invaginated cell shapes, but only actively contracting cable networks lead to the circular 

o : arc mo,phology ,ypical ** strongly adheri,,g cc " 8 ' In ,his case ' shape ™ d force distribllti0 " *" 

determined by local rather than global determinants and thus are suited to endow the cell with a 
robust sense of its environment. We also discuss non- linear and adaptive linker mechanics as well 

•i-H . 

^ ' as the relation to tissue shape. 

PACS numbers: 87.16.Ln, 87.10.-c, 87.17.Rt 



*I.Bischofs@zmbh.uni- hcidclbcrg.de 



i Ulrich.Schwarz@bioquant.uni-heidelberg.de 



1 



I. INTRODUCTION 



During the last decade, it has been increasingly realized that adherent cells actively ex- 
plore and respond to the geometry and stiffness of their adhesive environment, with dramatic 
consequences for fundamental cellular processes such as survival, proliferation, differentiation 
and migration [l, 2]. Using chemical inhibitors for both myosin motors and actin filaments, 
it has been shown that active contractility of the actin cytoskeleton is an indispensable re- 
quirement for the observed cellular response to the physical properties of the environment. 
Although more and more details are revealed regarding the molecular mechanisms under- 
lying these sensing processes 3ru>|, what is still missing is a systems- level understanding of 
how cells integrate the way force is generated, distributed, and sensed over the whole cell. 
Therefore theoretical models are required which describe how force is propagated through 
the actin cytoskeleton as a function of environmental geometry and stiffness. Because the 
actin cytoskeleton is an integral part of the sensing capabilities of cells, one expects that its 
mechanical properties have evolved to support these important cellular functions. 

Although it is a standard procedure in experiments to block force generation and propa- 
gation with chemical inhibitors or by RNA interference, it is very difficult to measure how 
force is distributed inside cells and between cells and their environment. Different experi- 
mental approaches have been developed to meet this challen ge. Traction force microscopy 



allows to measure the forces transmitted to the substrate 



10(1 . Recently this technique 



has been extended in such a way that also cell-cell forces can be estimated from cell-matrix 



forces 
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121 ] . However, it is important to note that many forces balance inside the cell and 



are not transmitted to the substrate, so the forces existing inside cells might be much higher 
than appreciated from traction force microscopy [l^|. Laser cutting allows to estimate forces 
from the mechanical relaxation after cutting load-carrying elements like microtubules in the 
mitotic spindle jl^, Ijj] or stress fibers in the actin c yto skeleton 
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19|. Laser ablation can 



be used for subcellular analysis of cortical tension 20]. However, these experiments only 
probe local relaxation events of prominent cytoskeletal structures and therefore might miss 
the global effects of more distributed and less visible structures. Micromanipulation can be 
used to distort the mechanical balance of the cell globally 2lM24j , but the resulting changes 



in force distribution can only be estimated indirectly from its effects, e.g. growth of focal 



2 



adhesions. To achieve a more systematic understanding, these experimental approaches have 
to be complemented by theoretical approaches. 

Because the mechanical properties of cells, extracellular matrix and tissue are strongly de- 
termined by filamentous networks of proteins like actin, tubulin, lamin, spectrin or collagen, 



mechanical networks are widely used theoretical models for cell and tissue mechanics 25l.l26|. 



One of the best studied cases is the red blood cell, whose shape and mechanics has been 



studied with network approaches in very large detail [27H35|. Modern computer power per- 
mits to simulate each of the roughly 10 5 spectrin links separately and with molecular detail, 
for exam ple using the appropriately parametrized force-extension curve of a semiflexible 
polymer [33j. In the limit of small extensions, these models usually reduce to Hookean 
networks. 

For adherent cells, the main structural determinant is the actin cytoskeleton, whose 
mechanics differs in several important aspects from the one of the spectrin network of red 
blood cells. In general, the molecular structure of the actin network is much less defined. 
Its most prominent feature in adhesion is strong contractility due to activity of myosin II 
motors. This observation implies that the mechanical links between the nodes of the network 
cannot be simple actin filaments, but have to be bundles of actin filaments crosslinked and 
tensed by myosin II motors. The simplest model for prestress in a mechanical network is 
the introduction of a finite resting length which is smaller than the typical extension of each 
link. Indeed, one- and two-dimensional spring networks with prestress are widely used for 
modeling cell migration 36l-l39|. 

Network models are conceptually very appealing due to their multi-scale nature: by 
changing the microscopic rules for the mechanics of the links, one can explore how the 
macroscopic behavior of the whole network changes. In particular, important biolo gica l 
effects like viscoelasticity of the links or coupling to diffusion fields can be incorporated 
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39l |. Spring networks offer the additional advantage that homogenization techniques can 



be used to arrive at continuum models. Recently, the interplay between force generation 
and the geometrical and adhesive properties of the environment have been addressed using 



the powerful framework of finite element models (FEM) 40|, |41] , which can be considered 
as the continuum limit of appropriate network models. Most FEM-models use constitutive 
equations which correspond to Hookean networks. 
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Although conceptually very appealing, modeling cell mechanics with Hookean networks 
does not reflect the fact that the actin cytoskeleton does not provide much resistance to 
compression. This is especially true for two-dimensional networks for cell adhesion and 
migration, because in this case the network might contract laterally, while the cytosol flows 
into the third dimension. In this situation, the network links do not behave as springs, but 
rather as cables, which are characterized by an asymmetric force-extension relation. There 
are several microscopic reasons for this effective behavior: not only do thin actin bundles 
easily buckle under load, they also tend to telescope in due to filament sliding and even 
to depolymerize once tensile stress is relieved. Cable networks have been successfully used 
to model the prestress-dependent mechanical response of adherent cells to local mechanical 



perturbations 



42] . The same model has also been used to describe how mechanical stress is 



propagated from the nuclear region through the cytoskeleton towards focal adhesions, where 
changes in load lead to changes in adhesion size 22]. 

One striking feature of strongly adhering cells is the fact that retracted contours often 
take the shape of circular arcs 43 -46j]. Although for cable networks the resulting shapes 
are strongly invaginated, it has been shown that the circular invaginations observed for cells 
pinned at discrete sites of adhesions can only be explained if an additional contractile force 
is introduced for each mechanical link j46|. This additional force in an actively contracting 



cable network does not vanish at the resting length and represents the fact that contractility 
arises mainly from myosin II motors, which in steady state operate close to a non- vanishing 
stall force. 

In this paper, we systematically compare the different network models introduced before 
for adherent cells (spring models, cable networks, and actively contracting networks) in 
regard to the predicted shapes and force distributions. Our main conclusion is that actively 
contracting cable networks share many crucial features with adherent cells. Due to their 
linear nature, actively contracting spring networks are equivalent to passive spring networks 
with a reduced resting length. Passive networks (both from springs and cables) have a well- 
defined reference state even in the absence of adhesion constraints and in general give similar 
results regarding shape and force distribution, which is determined mainly by global inputs 
like the spatial distribution of the adhesion points. In contrast, actively contracting cable 
networks do not have a well-defined reference state because without adhesion constraints, 
they contract onto a point. In this case, we find that shape and force distributions are 
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determined mainly by the local distribution of adhesion sites. The internal force distribution 
is constant in the bulk and strongly localizes to the contour, where forces jump by orders 
of magnitude. This motivates a detailed study of two contour models, which allow us to 
derive analytical predictions which we then compare with the results from the computer 
simulations. We also discuss how actively contracting cable networks can be extended to 
model non-linear or adaptive linker mechanics, and comment on the relation of our network 
models to tissue mechanics. 

II. NETWORK MODELS 
A. Link Mechanics 

Tissue cells adhering to discrete adhesion points on a flat substrate usually become very 
flat and therefore effectively two-dimensional. Only the nucleus, which rises in the middle, 
makes them fully three-dimensional, see sketches in Fig. [T^,,b. However, here we focus on 
lateral contraction and contour effects and thus the nucleus is expected to play a minor role. 
In the following we therefore restrict ourselves to two dimensions and model the cytoskeletal 
network as a two-dimensional mesh of mechanical links joined at N discrete nodes. Nodes 
are labeled with indices, e.g. % and j. Microscopically links and nodes may represent filament 
bundles and local accumulation of cross-linkers, respectively, but in a more general sense, 
these mechanical elements are simple representatives of an unknown network architecture 
which we model in a statistical sense. The network is subject to internal forces originating 
from molecular motor activity, F active , and the mechanical resistance of filaments to strain, 
Fmech, compare Fig. [Tb. 

We introduce three fundamentally different kinds of mechanical models for the network 
links. The simplest case is a Hookean spring network (HSN) composed of links with resting 
length Lq, which represent linear springs with spring constant EA/Lq, where E is the 
Young's modulus of the link and A its cross-section. The restoring force acting on a node i 
due to elastic strain in the link ij then reads: 

Fij,mech E AUijC^j , (1) 

where = (Rj — Ri)/Lij is the dimensionless unit vector along the link ij. Here, Ri and 
Rj specify the node positions, Lij is the length of the link and Uij = (Lij — L )/L is the 
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FIG. 1: (Color online) Sketch of the system, (a) Side-view of an adherent tissue cell, (b) Top-view. 
The cell is assumed to be adherent at four discrete dots. Its contour shows an invaginated shape 
resulting from a balance between bulk forces directed towards the cell center and boundary forces 
directed along the cell periphery, (c) The actin cytoskeleton is tensed by myosin II minifilaments, 
which actively contract the network with forces F ac ti ve . If the network links are strained, restoring 
forces F mec h appear, (d) Force-extension curve Fij(Lij) of a link ij in a simple Hookean network, 
a passive cable network, and an active cable network (in order of increasing dash lengths), (k = 
spring constant, T = motor stall force per length, Lq = initial link length, L c = critical length) 

strain in the link. The linear force-extension curve of a single link in a HSN is shown in 
Fig.dH dashed line with short dashes. 

The mechanical properties of the cytoskeleton are attributed mainly to the actin part. 
Actin is a semi-flexible filament prone to buckling under compression and thus behaves like a 
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cable, which can be stretched but not compressed. The Hookean assumption of a symmetric 
elastic response is therefore not valid. The mechanical properties of actin networks on a 
coarse-grained scale are more accurately described by assuming a finite resistance of filaments 
to tension, but no resistance to compression. The mechanical restoring forces originating 
from a link connecting two nodes % and j in the passive cable network (PCN) are therefore 
given by: 

Fij,mech = EAUij^ij, Lq < Lij (2) 
Fij,mech = 0, Lij < Lq. (3) 

We show the asymmetric force-extension relation of the PCN links in Fig. [TH as a dashed 
line with long dashes. 

Let us assume that molecular motors are homogeneously distributed in the network. 
Because they are arranged in a parallel fashion, their individual forces add up. We therefore 
assume that a link contracts with a force TL proportional to its initial length, where T is 
force density per length. For computational simplicity, here we assume that this force does 
not change as the filament contracts, although in practise, it might well be that the line 
density rather than the total number of active motors is constant. For a link ij we therefore 
have: 

Fij^active T Lq€^ , (4) 

where T > is the tensile force per initial length applied by the motors. The finite force at 
zero length is unphysical and we avoid it by introducing an additional rule such that force 
is diminished if two neighboring nodes come closer to each other than some small distance 
L c < L : 

Fij,active TLq — Gjj, Lij <C L c . (5) 

Lie 

Combining PCN and active contraction, we obtain what we call the active cable network 
(ACN), compare the solid line in Fig. [TH: 

Fij = (TL + EAuij)eij, L < Ly, (6) 
Fij = TLoCij, L c < < Lq, (7) 
Fij = TLQ—^-eij, L^ < L c . (8) 

Neither the detailed choice of L c nor the assumption of a linear force reduction below L c are 
crucial for our results. 
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Active contraction can also be combined with the HSN. However, this simply shifts the 
straight dashed line in Fig. [T]i, i.e. this reduces the resting length L to 

4 = A, (i - 1?) (9) 

as long as TL < EA. Therefore active contraction does not change the basic definition of 
the HSN and therefore in the following the actively contracting HSN is not discussed further. 



B. Mechanical Equilibrium 

For nodes within the network, the total force exerted on a node i is the sum of all forces 
applied by neighboring nodes j: 

£ = £4- (io) 

3 

In mechanical equilibrium, the force on each non-adherent node has to vanish, Fi = 0. The 
existence of adhesion sites is modeled by fixing the positions of the adherent nodes. Thus 
the adhesion site geometry will enter through the boundary conditions. 

To reduce the number of parameters, we scale all lengths with respect to L , e.g. £ = L/L . 
All forces are scaled as / = F/EA. We define the ratio of active to elastic forces as 

^ = — • (11) 
EA y ' 

For an ACN we therefore may rewrite the forces acting on a node in the non-dimensionalized 
form as: 

fij = {uij + T)eij l<£ij, (12) 

fij = T&ij — ^ij — 1) (13) 

t% ■ _ 

fij = T ~ ei i ^ (1^) 

In the computer simulations we use i c = 10 -3 . 

Mechanical equilibrium requires the forces on each non-adherent node to vanish 

fij = V non-adherent nodes i (15) 

i 

with the summation j over all neighbors. For a two-dimensional network of N' non-adherent 
nodes, the system of equations (fl5]l consists of 2^' coupled non-linear equations. If the lhs 
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of system (Tl5|) consisted of arbitrary functions of the r*j, the method of choice to solve it 



would be the Newton-Raphson method 
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However, since the lhs of system f JT5|) is a force 



which has a potential, it is also a 2A r '-dimensional gradient vector. Therefore we solve the 
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We stop 



minimization problem for the potential with the conjugated gradient method 
iterating as soon as the force on every node (except the periphery nodes) is smaller by at 
least two orders of magnitude than the smallest link force. 

C. Parametrization 

Due to its dynamic and multiscale organization and the known limitations of microscopy, 
a detailed model of the actin cytoskeleton is currently out of bounds. In the face of these 
uncertainties, our model is not meant to represent the details of the organization of the 
actin cytoskeleton. Nevertheless, for practical purposes it is helpful to parametrize our 
model using some benchmark values for the actin cytoskeleton. 

The elastic modulus of an actin filament, which has cross-section area Afn = 18.8 nm 2 , 
was experimentally found to be Efu = 2.8 GPa |49J, while typical values for stress fibers 
are a radius around 100 nm (corresponding to A jib = 31416 nm 2 ) and an effective modulus 



of E fib = 1.45 MP a 



50| | . Hence, the Young's modulus of stress fibers is three orders of 
magnitude smaller than the one of single actin filaments. This suggests that cross-linkers 
like a-actinin and myosin II are the main contributors to elasticity and not the actin filaments 
themselves. However, the values for the one-dimensional modulus, EfuAfu = 52.6 nN and 
EfibAfib = 45.6 nN, are effectively very similar, so the one- dimensional modulus is expected 
to be of the order of 50 nN. 

The mesh size of the cytoskeleton is expected to be typically around L Q = 100 nm. This 
is an intermediate value introduced i n |42l| based on experimental observations of the actin 



cytoskeleton in adherent endothelial [51] and fibroblast cells 52J. For the active force, we 



estimate that around 1.000 myosin II motors are active in one effective link. With a stall 



force of 2 pN per motor head 53j, we have T = 2 • 10~ 2 nN/nm for the motor force per 
length. Using Eq. (ITT]) and EA = 50 nN, we estimate r = 0.04 for the active tension in the 
network. 
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(a) K d = ioe =>| (b) d = io* H 



FIG. 2: Tension- free reference states, (a) Square network with link length I = 1 and side length 
<i = 101. (b) Triangular network with the same t and d. 

III. RESULTS 

A. Equilibrium Shapes 

In the following we analyze how the different types of networks act under tension for given 
adhesion constraints. We first discuss passive networks. In order to understand the role of 
geometry, as tension-free reference states we use both a square and a triangle as shown 
in Fig. [2j If the corners are fixed and the network is set under tension by reducing rest 
length, invaginated shapes appear as shown in Fig. [3] for Hookean spring networks (HSN) 
and passive cable networks (PCN). Contraction is quantified by m = {t — £o)/£. In Fig. [3l 
the results for HSN and PCN are the same for the square shape, Fig. [3^l,c, because these 
two kinds of networks behave identical as long as all links are tensed. However, the results 
are different for the triangle shape, Fig. [3jD,d. In this case, the PCN gives a significantly 
flatter contour due to the missing response to compression in the thin extensions leading to 
the adhesion points. The most prominent examples for compressed links in Fig. [3p,d are 
indicated by arrows. 

In order to understand our numerical results in more detail, we first note that the HSN 
with triangular network topology has a well-defined continuum limit, in which it corresponds 
to a two-dimensional sheet with isotropic linear elastictiy 0|- The ^ wo correspond- 
ing elastic constant are a Young's modulus of 2/c/y / 3 and a Poisson's ratio of 1/3, where 
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FIG. 3: (Color online) Tensed Hookean (HSN) and passive cable networks (PCN) with th = 0.2. 
The colorbar gives the dimensionless force /. (a) Contraction of a HSN in the square reference 
state, (b) Equilibrium shape of a HSN in triangular reference state, (c) PCN in square geometry, 
(d) PCN in triangular geometry. For the triangular reference shape, HSN (b) and PCN (d) differ 
in the amount of compression in the thin extensions (marked by arrows). 



k = EA/Lq is the spring constant of the links. The HSN with simple cubic topology does 
not have such a rigorous limit, but in our context it works in a similar way as the triangular 
lattice. Therefore similar results as obtained here for the HSN are also obtained with con- 



tinuum elasticity theory applied to two-dimensional cell shapes |40|. Without any adhesion 
constraint, the HSN contracts isotropically to a finite size, i.e. the network is uniformly 
scaled and has a new side length d! = 10£q. This shape we call the unconstrained reference 
shape and it is key to understand the results for HSN. The same shape as shown in Fig. [3^ 
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FIG. 4: (Color online) (a) Boundary line of the HSN or PCN for the square shape for network 
tension th = 0.2 (lower set) and th = 0.6 (upper set). The straight line shows the unconstrained 
reference shape, while the three curved lines represent different initial link lengths, namely I = 
1, 0.1, 0.02 (from bottom to top within one set), (b) Boundary line of the active cable network 
(ACN) for the square shape for r = 10~ 3 (lower set) and r = 10 -2 (upper set). In this case, 
no unconstrained reference shape exists. The three lines again represent the initial link lengths 
£ = 1, 0.1, 0.02 (from top to bottom). 

results if the network contracts away from its initial state under adhesion constraints or if 
the network starts from its unconstrained reference state and its corners are dragged to the 
desired adhesion points. This explains the main feature of the force distribution shown by 
the color coding in Fig. [3^, namely the strong localization of stresses and strains to the 
regions around the adhesion points. With the amount of tension used here, the network 
can attain its unconstrained reference state away from the adhesion points and therefore its 
contour is essentially flat in the middle parts. 



12 



In Fig. [4^ we directly compare the calculated network shapes for the HSN for the square 
geometry to the unconstrained reference shape. In addition, we demonstrate the role of the 
link length i. As explained above, in this case the PCN gives the same results. For the small 
value of tension, th = 0.2 (lower set), the contracting network can reach the unconstrained 
reference shape over a large region where it is therefore essentially flat. The smaller £, the 
faster this contour is reached due to an increased force density along the contour. For the 
large value of tension, m = 0.6 (upper set), the unconstrained reference shape is not reached 
by the contracted network and it stays non-flat along its whole contour even for rather small 
values of I. 

Next we turn to actively contracting networks. Figs. |5h-d show the equilibrium shapes of 
an active cable network (ACN) with reference state from Fig. [2k. and increasing tension r. As 
tension increases, the shape becomes more and more invaginated, until it collapses onto the 
zero-area network in Fig. [5H. This network basically consists of a centrally contracted region 
which is connected to the adhesion points by long extensions. We therefore call this network 
the center tree (CT). Note that this network still retains aspects of the two-dimensional 
network, because an effectively one-dimensional structure would collapse onto the so-called 
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Steiner tree of minimal length, which for a square shape is not four-fold symmetric 

In contrast to the passive networks, where tension CCR.S6S cLS the unconstrained reference 
state is reached, for the ACN no such unconstrained reference shape exists and without 
adhesion constraints the shape would collapse onto a single point. This explains why flat 
parts are not observed in the contours of the networks shown in Fig. [5^-d. This is also 
demonstrated in Fig. 0b, where we show contour shapes for r = 10~ 3 (lower set) and 
t = 1CT 2 (upper set). Like for the HSN in Fig. [4^, increasing lattice constant i increases 
the force density along the contour and therefore leads to stronger invagination. However, 
in marked contrast to the HSN, no flat parts appear in the contour as no unconstrained 
reference state exists. 

Fig. E^-d also shows that the formation of inward directed arcs now corresponds to a 
much more inhomogeneous density distribution of filaments: in the bulk of the network 
the distance between nodes and thus the filament density remains unchanged, while at the 
edges filaments start to bundle strongly along the edge. The color code in Fig. [5] shows that 
stress is strongly localized at the periphery. In the interior, the only forces acting are the 
motor forces r which balance each other at every node. At the periphery, the force jumps 
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0.001 0.01 0.1 1 

FIG. 5: (Color online) Contraction of an active cable network (ACN) with square shape, (a) and 
(b) Contraction of the network leads to arc formation, (c) With increasing tension tubes form near 
the adhesion points. Note the bundling of filaments at the edge, (d) For very large r, the network 
collapses onto the center tree, (e) The same situation as in (a) but with t = 0.5. The contour is 
equal to that in (a) and the boundary forces are the same, (f) Square Voronoi network with 212 
nodes and 316 links. Note the regular contour and the strong stress localization in the periphery. 
Tension values are (from (a) to (f)): r = 10~ 3 , 10~ 2 , 10" \ 1, 1CT 3 , 1(T 3 . 
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up from r to much higher values r + EAA£j/£, compare Fig. [5^. The forces are largest close 
to the adhesion points and decrease towards the center of the boundary. For large tension 
(r > 10~ 2 ), tubes are formed near the adhesion points and the stress distribution along the 
contour becomes more inhomogeneous. 

Fig. shows the effect of changed discretization for the same tension value as Fig. 
for the complete network. Fig. [5]F demonstrates that for ACN, shape and force values do 
not depend significantly on network topology. As an instructive example here we use a 
disordered network topology obtained by a Voronoi construction. Even the presence of 
relatively large elements in the discretization does not change the invaginated shape feature 
of the contracted network. We conclude from Fig. and Fig. O that ACN are surprisingly 
robust in regard to the details of the network topology. In general, similar results as obtained 
in Fig. [5] for the square shape are also obtained for other initial cell shapes. 

In Fig. [61 we investigate another striking property of ACN, namely its robustness in 
regard to addition of new adhesion points. In Fig. |6K,b we show the equilibrium shapes 
of the HSN from Fig. |3k and the ACN from Fig. |5b with one adhesion point added in the 
middle of the bottom line. In the case of the ACN model, this change at the bottom of the 
network has little influence on the positions of nodes not directly connected to the bottom 
line. In contrast, for the HSN the additional adhesion point affects the shape of the opposite 
arc, becoming more curved in the center. In Fig. [6fc we plot by which distance S the nodes 
in the vertical middle line are pulled down in the negative y-direction upon addition of the 
new adhesion point. The plot of log 5 versus y is not smooth for numerical reasons, but 
clearly shows that the effect decays much more rapidly for the ACN versus the HSN. In 
addition to square and triangular network topologies, here we also show results for rotated 
square (diamond) and hexagonal networks. Intriguingly, stress in the contour behaves very 
differently, compare the color coding of Fig. (6^,b. While in the HSN stress in the bottom 
line stays approximately the same, in the ACN it decreases to half its value, indicating a 
strong effect on contour forces. 

In summary, ACNs behave very differently from HSNs (and therefore also from the mostly 
equivalent PCNs). Roughly speaking, they respond more locally than globally. They are 
more robust in regard to network topology and adhesion geometry and show strong local- 
ization of the stress to the periphery. Addition of new adhesion points leads to little global 
change, but to a strong change in local stress distribution. Because for ACN shape and 
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FIG. 6: (Color online) Contraction of networks with square shape with an additional adhesion 
point in the middle of the bottom line, (a) HSN or PCN model with th = 0.2. (b) ACN with 
r = 10 -2 . (c) Relative displacement 8 of nodes with x = 5 (vertical middle line) with and without 
the additional adhesion point at the bottom, y is the node's y-coordinate in the initial network. 
The four top lines correspond to the HSN while the four bottom lines represent the ACN. Different 
symbols show different topologies (square=D, diamond=0, triangular=A, hexagonal=O0 

contour stress seem to be mainly determined by the local adhesion geometry, they will now 
be analyzed in more detail. 

B. Contour shape and tension-elasticity model 

In contrast to HSNs, the contour of ACNs appears to be more circular. Indeed, a circular 
arc morphology has been noted before for the shapes of cells adhering to micropatterned 
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FIG. 7: (Color online) Arc fits, (a) Arc analysis for a tensed HSN or PCN. From bottom to top: 
TH = 0.1 (+), 0.2 (*), 0.4 (x), 0.6 (O), 0.8 (□), 1 (0). Symbols: Contour. Lines: Straight fits of 
the linear contour parts (the three bottom lines), least square fits of round contour parts to arcs 
with constant curvature (the three top lines). * corresponds to the bottom line from the networks 
shown in Figs. [3^ and c. (b) Contour analysis for the bottom line of an ACN with r = 10~ 3 (x), 
10~ 2 (+), 2.5 -10" 2 (*), lO" 1 (O), 2.5 -10- 1 (□), 5-10- 1 (A) from bottom to top. Symbols denote 
node positions of the network arcs while the lines are least square fits of the contours to arcs of 
constant curvature, x, +, Q> an d dashed line correspond to the bottom lines from Figs. [5^i-d. 
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substrates and therefore this shape feature is an important motivation to study ACNs 46 ]. 
We now investigate this important aspect in more detail. In Fig. [7^ and b we show contours 
of the HSN and ACN from Fig. [3] and Fig. \5[ respectively, together with fits to straight 
lines or circular arcs. As we increase tension, we observe more invaginated shapes, compare 
Fig. HJ For the HSN in Fig. [7^,, small tension allows the network to reach the unconstrained 
reference shape and therefore the best fit to the middle part of the contour is a straight 
line. For larger tension, the unconstrained reference shape cannot be reached anymore and 
circular shapes become better fits. This crossover is in marked contrast to the ACN from 
Fig. [7b, where circular arcs fit very well for all values of r. For large r, the overall contour 
starts to deviate from the perfect arc shape because the networks starts to collapse into 
tubes near the adhesion points. However, locally (in between the tubes) the contour stays 
circular. Another difference between the two network types lies in the observation that for 
HSN, network shape strongly depends on lattice constant £, while for ACN, the equilibrium 
contour is relatively independent of £. 

It has been argued before that the circular arc shape feature of the ACNs can be explained 
by an analytical theory, the tension-elasticity model (TEM) 46[j. For clarity, here we repeat 
this analysis and compare it in detail with our network simulations. Because ACNs do not 
propagate compression and the motor forces represent a constant pull in the network, in 
the TEM the bulk contractility is modeled by a structure-less surface tension a. However, 
elasticity is crucial to understand how the contour reacts to the internal pull. Therefore the 
elastic nature of the mechanical network is represented by an elastic line tension A, which 
prevents the contour from collapsing under the inward pull of the bulk network. This line 
tension is written as 

A = EA^^ (16) 

where L is the contour length and L is its resting length. Note that we use the same 
value EA like for the single links because the elastic line tension will be dominated by the 
contribution from the most peripheral line of links. We further assume Lq = ad, where 
d is the initial {spanning) distance d between two neighboring adhesion points (which in 
the simulations above has been chosen to be 10) and a is a dimensionless resting length 
parameter (compare Fig. [8^l for a schematics). In the following we restrict ourselves to 
a = 1, that is we assume that a completely relaxed contour is straight, but without internal 
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FIG. 8: Schematic representation of the two contour models, (a) In the tension-elasticity model 
(TEM), an isotropic surface tension a pulls the contour in along the normal direction, while the 
counteracting line tension A acts along the tangential direction, r is arc radius, L is contour length, 
and d is spanning distance, (b) In the elastic catenary model, the inward pull is vertical and thus 
leads to an inhomogeneous line density of force along the elastic contour, eft denotes the tangential 
angle. 



tension. This implies that we neglect the contribution of the active contractility in the 
periphery to the line tension. 

The relation between surface tension a and network tension r, which depends on network 
topology and discretization, can be obtained numerically. For this purpose, we simulate the 
pulling of a rectangular sheet of network. Surface tension a then follows as total force on 
the pulling boundary divided by its width. Due to this normalization, the result does not 
depend much on link length £, thus we use £ = 1 for the simulations. For all considered 
network geometries, we find a linear relation: 

^square 0.9907 ■ T~ S q Uare} (1^) 

^triangular 1.6892 • Ttriangular i (l^) 

^diamond 1.0867 • T&iamondi (19) 

^hexagonal 0.5517 • T}i exa g 0na i. (20) 

The constant for the square lattice is close to 1 because here all links pull essentially per- 
pendicular to the boundary. 
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FIG. 9: (Color online) Relation of arc radii to internal tension of ACN for different dot distances 
and comparison to the tension-elasticity model (TEM). Symbols denote simulation results (□ 
corresponds to square, to diamond, A to triangular, and O to hexagonal topology), the solid 
line is the numerical solution of Eq. (|'24|) and the dashed line the analytical result Eq. (|'25|) . Side 
lengths are (from bottom to top): d = 10, 19, 31, with critical tensions a c « 0.114, 0.060, 0.036. 

Given the forces assumed by the TEM, one can derive the shape of the contour from the 
force balance. While the surface tension a acts in the direction of the normal n, resulting 
in a pulling force F = an, the line tension A acts in the tangential direction t (compare 
Fig. [Hk). Because the elastic line tension is a global quantity, it does not vary with the 
contour length s and therefore the contour tension is T(s) = Xt(s). Then the force balance 



the radius of curvature. We thus conclude that the TEM predicts circular arcs with a radius 




n/r with r being 



(21) 



A 



(22) 



r = — 
a 
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Although this results looks like a simple Laplace law in two dimensions, it is more compli- 
cated, because the arc radius r will depend on global properties like spanning distance d 
through the elastic line tension A from Eq. (j!6p . 

In order to arrive at an expression for arc radius r as a function of adhesion geometry 
and network tension, we use the trigonometric relation 



to replace contour length L by spanning distance d, compare Fig. [SJa,. In combination with 
Eq. (|T6|) (in dimensionless form) and Eq. (122]) . this gives 



Since this equation cannot be solved analytically for r, it has to be solved numerically 
for given values of d and a. For geometrical reasons, r must always be larger than d/2. 
Therefore, a critical a c exists above which Eq. (1241) cannot be solved anymore. For small 
values of a, the invagination is small and one can expand the geometrical relation md/r <C 1. 
This leads to the analytical result 



In Fig. [9] we compare the results from computer simulations for r over a large range of 
network tension r to the results of the TEM with the corresponding range of surface tension 
a, both for the numerical solution of Eq. (|24|) and the analytical solution Eq. (125]) . We note 
that the predicted power law behavior applies over a very large range of tensions, and only 
breaks down at very large tension a > 10 -2 , compare Fig. [9j The inverse relation between r 
and u represents a modified Laplace law for ACNs and thus demonstrates that the concept 
of an isotropic surface tension works well to explain cell shape. With the linear relation 
between a and r, this implies that r ~ r~^. 

We also find excellent agreement between computer simulations and TEM upon variation 
of spanning distance, compare Fig. [9j Thus the elastic effects mediated by the spanning 
distance d are captured well by the concept of an elastic line tension. In summary, the 
analytical TEM results in a surprisingly good description of the contour shape of ACNs. As 
we will discuss in the next section, however, agreement is less good regarding contour forces. 




(23) 



r = — 




(24) 
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r = 24 3rf3<7 3. 
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FIG. 10: (Color online) Force distribution in adherent networks, (a) Force of vertical links which 
cross the straight line y = 5.5 in the HSN from Fig. [3^ (top) and the ACNs from Fig.[5h.,b (bottom, 
middle), (b) Force in the bottom line links of an ACN with a = 10~ 3 for different lattice constants 
(symbols). From top to bottom: £ = 2, 1, 0.5, 0.2. The straight line gives the line tension obtained 
via the TEM, the curved line follows from Eq. (|3ip . In both, (a) and (b), on the x-axis we have 
the x-coordinate of the center of mass of the links, (c) Forces on adhesion dots exerted by the 
ACN (top) and HSN (bottom). Simulation results are shown as dots, while lines give the power 
law fits. For the ACN we obtain / = ar b with a ~ 6.16 and b ~ 0.783 and for the HSN / = ar\j 
with a ~ 4.51 and b ~ 1.26. (d) Power law fits of adhesion dot force vs. surface tension. ACN and 
elastic catenary results (top line, collapsed) and TEM (low line). 



C. Contour forces and elastic catenary model 

We now discuss the forces resulting from our computer simulations. In the simulated 
networks of Fig. |3] for HSN and Fig. |5] for ACN, the stresses in the network are coded by 
colour. In Fig. [TUh we plot these stresses along a line horizontally crossing the networks 
with square shapes slightly above the middle line. In the HSN, stress gradually decays into 
the sample, while for ACN, it jumps up at the periphery. Fig. [TUh also shows that forces in 
HSNs are much larger than those in ACNs with a comparable equilibrium shape. The stress 
distribution in the boundary of an ACN depends on the lattice constant I of the network, as 
shown in Fig. [TUb . The contour forces are minimal in the middle and increase towards the 
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sides. At an adhesion site, all network forces add up to the overall force being transmitted 
to the substrate. For both passive and active networks, the force which is exerted on an 
adhesion site follows a power law as tjj (r) is increased, see Fig. [TUb . For r = th, adhesion 
force is much smaller in the HSN than in the ACN. For active and passive networks of a 
comparable shape, however, e.g. th = 0.2 and r = 0.01, we observe the opposite behavior. 

Although on an absolute scale the variation is not very strong, Fig. [5b and Fig. [TUb 
both demonstrate that for ACNs under large network tension, peripheral force varies along 
the contour. In contrast, the tension-elasticity model (TEM), which is very successful in 
explaining shape, predicts homogeneous force A = ra along the boundary. Fig. \5jp suggests 
one reason which could explain this discrepancy. For ACN, the links essentially telescope in 
under contraction and therefore their density along the contour varies for strong curvature 
along the contour. This suggests that in order to explain the spatially varying force in 
the contour, one has to revisit the assumption of an isotropic surface tension a creating a 
homogeneous force density along the contour. 

As an alternative to the TEM, we now investigate another analytical model, which in- 
corporates the effect of varying link density, namely the elastic catenary 55[ . In the elastic 
catenary, the pulling force on the elastic contour is not along the normal, but along the 
direction perpendicular to the original contour, similar to the situation in networks with 
square shape and square topology, compare Fig. [8b. Due to the linear elasticity in the con- 
tour, the line density of links along the contour varies in proportion to the contour tension. 
In dimensionless units, a length element is expanded to length 1 + T(s) in the presence of 
tangential contour tension T(s) = T(s)t(s). With an initial inward force per unit length a, 
the effective force density is therefore a/ (1 + T(s)). In contrast to the TEM, we now assume 
not only a heterogeneous force distribution, but also a vertical pulling direction. Therefore 
the pulling force is F(s) = (0, —<r/(l + T(s))). The force balance again reads 

f + (26) 

but now the solution is more difficult than for the TEM. Because the tangent is normalized, it 
can be written as t = (cos0(s), sin0(s)), where <j)(s) is the tangential angle along the contour 



56J, compare Fig. [8b- Different from the TEM, we now have to solve two equations: 



^(T( S )cos(0( S ))) = (27) 
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>W s «) = TTTw- (28) 

Eq. ([27)1 can directly be integrated, leading to Tcos(0) = const = A c , while Eq. (128|) can 
then be solved via the substitution tan(0) = sinh(p). This gives 

A A 2 

x(p) = —p H — - sinh(p) + x , (29) 
a a 

A A 2 

Vip) = — cosh(p) + cosh 2 (p) + y . (30) 

The integration constants Xq and i/q are determined by the positions of the adhesion sites. 

While in the TEM we assume the line tension A to be constant along the whole boundary 
line, for the elastic catenary only the x-component of tension, A c , is constant. Its value 
can be determined numerically from the above equations for given d and a. One can show 
that for small surface tension a, the contour becomes parabolic with a radius of curvature 
r = A c (l + A c )/cr. For cr = 10~ 2 the relative deviation to the prediction of the TEM is only 
4%, that is in this regime, the elastic catenary model leads essentially to the same result 
as the TEM with circular arcs. However, in contrast to the TEM, this model predicts a 
spatially varying boundary tension of 

T{p) = A c cosh(p) . (31) 

The curve without symbols in Fig. [TOb shows that this model qualitatively predicts the 
observed minimum in the stress distribution. The force acting on an adhesion dot is predicted 
to be 

/ = V2Ae(l + Binh(po)). (32) 
Fig. [TOU shows that this prediction is quite accurate. 



D. Strain Stiffening 



So far, we have treated the network links as cables or springs, in which force increases 
linearly with elongation. This implies that the elastic modulus of the network is indepen- 
dent of strain. However, it is well known that the cytoskeleton of cells shows a strong 
increase of elastic modulus with strain {strain stiffening) 57j. This property of cytoskeletal 
networks has been reconstituted in vitro with crosslinked actin filaments of physiological 



length, which is of the order of their persistence length (around 1 /jm) 



Under these 
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FIG. 11: (Color online) Contracted networks with worm-like chain (WLC) mechanics, (a) PCN 
with t h = 0.47. (b) PCN- WLC with t h = 0.47 and l n = 0.5. (c) ACN-WLC with r = 1 and 
t n = 0.5. (d) Maximum invagination h of linear (upper curve) and non-linear (lower curve) ACN. 



conditions, strain stiffening results mainly from the mechanical properties of the crosslinkers. 
For larger filament lengths and stiff crosslinkers, there exists another mechanism for strain 
stiffening, namely the non-linear force-extension curve of seminexible polymers as described 
by the worm-like-chain (WLC) mod el |59l| . The WLC-model has been used before to model 
seminexible biopolymers like DNA [60(, actin 6l|, |62[ and spectrin 33|, |34]. Because our 
network model is especially suited to study the effect of link mechanics, we now study strain 
stiffening based on the WLC-model, although for the physiological relevant case, a model 
for crosslinker mechanics might be more relevant. 
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While the WLC proper has vanishing resting length, here we combine it with a finite 
resting length to also include the effect of compression. Thus we use the non-linear WLC 
model to describe the mechanical links as they are tensed away from their reference state, 
while the compressed state is modeled as above (linear response for springs and no response 
for cables). Complementing Eqs. (J2J) and for the PCN, we get in dimensionless form 



Force is now given as multiples of ksT '/ Lq, while length is again scaled with Lq. Strain is now 
defined as Uij = (£ij — l)/£ n , in which l n gives the dimensionless difference between maximal 
extension and reference length. Effectively there is only one difference to the original model, 
namely the additional term which diverges if the strain it^ approaches 1. Without strain, 
this term vanishes. For ACNs, the original model is extended in the same way. 

We choose l n = 0.5, that is the maximal extension is 1.5 Lq. For small tension, r# < 0.2, 
the non-linearity does not affect the shape of the PCN much and we observe the same 
invagination as in Fig. |3j Fig. [TTb.b shows for the PCN a comparison between linear and 
WLC- networks for a large value of tension, t# = 0.47. Obviously the strain-stiffened network 
shows a much larger resistance to invagination. We also note that forces are two orders of 
magnitude larger in the non-linear model. ACNs are affected less by the non-linearity, as 
shown in Fig. [TTc . Here we use r = 1 and again i n = 0.5. Comparison with Fig. [5H 
reveals that the ACN collapses to a lesser degree than without strain stiffening. This can 
be quantified by the arc height h(r), defined as the maximum distance between initial 
and current edge in the equilibrium shape, h is significantly reduced by the non-linearity, 
Fig. EH. Thus much higher motor forces are needed to reach the collapsed state. 

In Fig. [T2k . the contour of the strain stiffening PCN is analyzed in more detail. Circular 
and linear fits are shown as lines. The bottom two are given by straight lines. For t# < 0.3 
the contour of the PCN- WLC is qualitatively the same as that of the linear one, shown in 
Fig. [7b. At th = 0.3 the contour cannot be fitted well by circle or line. If th is increased 
beyond 0.3, the network does not contract any further, but again expands outward. This 
surprising effect does not occur for ACNs. For th > 0.4, the arcs appear to be circular. 
The ACN- WLC contour, Fig. [12b, only differs little from the linear ACN contour, Fig. [7b. 
Arcs are always circular (except at the regions where tubes form). With increasing r they 





(33) 



fij Eij ^ 1 . 
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FIG. 12: (Color online) Arc fits for HSN and ACN with non-linear links, (a) Contour of the 
HSN-WLC from Fig. HQ). Symbols belong to different values of t h : 0.1 (+), 0.2 (*), 0.3 (x), 0.43 
(□), 0.47 (0)- correspond to the bottom line of the network from Fig. Illb . (b) Contour of the 
ACN-WLC from Fig. [TTb . Symbols are bottom line node positions, while the lines are circular fits. 
Motor force values are r = 10~ 2 (x), 2 • 10~ 2 (+), lO^ 1 (*), 2.5 • KT 1 (O), 5 • KT 1 (□), 1 (A). 
Note, A gives the bottom line from Fig. [TTb . 

continuously move inward. Comparison with Fig. [T^i reveals that radii typically are larger 
in the non-linear case. 
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FIG. 13: (Color online) Equilibrium shapes for adaptive networks, (a) Square HSN/PCN network 
with ea ~ 9.40 for the peripheral links and ea ~ 1.75 for the innermost links, (b) Square ACN 
network with ea ~ 7.10 at the periphery and ea ~ 2.00 for the internal links. Parameters are: 
Hi = 10, f = 0.5 (in (a)), r = n = 10" 2 , [eo]i = 10, f = 0.1 (in (b)). 

E. Link Adaption 

Adherent cells are known to strongly adapt their cytoskeleton to the physical properties 
of their environment. During recent years, it has become clear that the actin cytoskeleton 
tends to reinforce under load. In addition, mechanical loading of adhesion contacts leads to 
regulatory signals which increase myosin motor activity inside the cell. Network models are 
especially suited to study these biologically important effects in a theoretical framework. In 
the following, we will investigate which changes occur in the network if the elastic constant 
EA and the force density T resulting from myosin II activity are increasing with load. 

For simplicity, we assume that both EA and T first increase with force in a linear fashion 
and then saturate at constant values, which is the simplest assumption for a process based 
on enzymatic regulation: 

EA(F) = [EA] + [EA]i 1 £ w , (35) 
T(F) = To + 71-^— , (36) 
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where the force scale Fq determines when half the maximal increase has been reached. We 
again use dimensionless parameters. Forces EA, T, F are measured in units [-EM]o, i.e. 
we define r = TL /[EA} (the same for r and Ti), ea = EA/[EA} (the same for [ea]i, 
[eajo = 1) and / = F/[EA] (the same for / ). 

Fig. [13] demonstrates that the effect of adaptation is fundamentally different for passive 
versus active networks. Here stiffer links are represented by thicker lines. Fig. [T3h shows 
the result for a HSN (as shown by Fig. [3] for the square shape a PCN gives the same 
results). In this case, only Eq. (135]) must be considered. While the peripheral links show 
the largest values of ea, the rigidity decays smoothly from there into the bulk. This is 
strikingly different for the ACN shown in Fig. [T5b . Here Eqs. ( 13 5 j) and (I36p have been 
used. We find that the stiffness is strongly localized to the periphery, as found before for 
the internal stress. Thus for passive networks the adaption response is spatially continuous, 
while for actively contracting networks, it is strongly localized to the rim. This nicely agrees 
with experimental observations that strong peripheral actin bundles typically line the cell 
contour 46] . In particular, our model suggests that this effect is strongly determined by the 
mechanical properties of the underlying networks. 



F. Relation to Tissue Shape 



Tissue contraction with discrete pinning sites is very similar to cell contraction since 
adherent tissues also show invaginated arcs 46j, [63j. However, in this case the spatial di- 



mensions of the pinning objects tend to be relatively large. In order to include this effect, 
we have simulated a circular network contracting around four circular dots of finite size, see 
Fig. [T4"a . All nodes on the circles are fixed in space. Motor force r is increased stepwise and 
nodes coming closer to each other than £ c = 0.01 are glued together and in the following 
act as one 22J. Relative dimensions are taken from 63J], where a microtissue tethered to 
four cylindrical posts is analyzed. With the ACN we are able to reproduce the typical arc 
morphology of the contracted microtissues, see Fig. [T3b. We note that this is not possible 



with a FEM approach, as this leads to flat contours as for HSN 
model both for cells and tissues. 



63l |. Thus ACN are a useful 
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(a) 



(b) 



FIG. 14: (Color online) Contraction of an ACN anchored to dots of finite size, (a) Initial situation. 
Round tissue with radius r% adherent to 4 round dots of finite radius which form a square with 
side length dd- (b) Contracted tissue. Parameters are: rt = 38, = 4, da = 36, r = 0.1. 

IV. CONCLUSIONS AND OUTLOOK 

Motivated by the network nature of the actin cytoskeleton and its effectively 2D organiza- 
tion in mature adhesion to flat substrates, we have modeled adherent tissue cell contraction 
by 2D network models. The main aim of this work is to achieve a detailed comparison of 
the shapes and force patterns for different network types, namely Hookean spring networks 
(HSNs), passive cable networks (PCNs), and active cable networks (ACNs). 

The shape of a HSN can be understood best by considering the shape of its unconstrained 
reference shape. If tension is not too large, the network contour follows the unconstrained 
reference shape at regions sufficiently far away from the adhesion sites. Closer to the adhesion 
sites, the network deforms and stress and strain accumulate. In contrast, the ACN does not 
have an unconstrained reference shape and without adhesion constraints would contract 
into a point. Therefore no signature of the unconstrained reference shape (like flat parts 
for a square-shaped lattice) appear in the contour. Because it does not resist compression, 
stress and strain are not propagated much into the network and are strongly localized to 
the contour. 

One of the most striking difference between the different network types revealed by our 
analysis is the fact that in passive networks, local changes to the adhesion geometry changes 
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the network globally. This is in marked contrast to the active network, where the addition 
of local adhesions has only a local effect on the boundary. However, in this case the change 
in spanning distance has a large effect on the stress in the contour, as predicted by the 
tension-elasticity model (TEM). The TEM is especially suited to quantitatively predict the 
shape of an ACN, namely the circular arcs observed between neighboring adhesion points 
and the scaling of their radius with spanning distance and surface tension. 

Despite this success, the TEM does not capture all aspects of the network model. While 
the TEM assumes constant contour tension, the computer simulations reveal that tension 
varies along the contour. An elastic catenary theory qualitatively predicts that tension 
decreases towards the middle of an invagination due to local changes in link density along 
the contour. However, it does neither predict the quantitative details of the contour stress 
nor the circular arc morphology. 

In the case of very large tension, both network types develop different features. For 
the passive networks, the invaginations tend to become more round, as the unconstrained 
reference shape becomes so small that the contour cannot reach it anymore. In contrast, 
the active network develops straight features, because the network collapses into tubes at 
the adhesion points. Indeed the formation of tubes has been observed experimentally and 
eventually leads to pearling through a Rayleigh-Plateau instability 44|. The region between 
the tubes always stays circular for ACN. 

Network models are ideally suited for multi-scale modeling because physical properties 
can be easily added on the level of single link and lead to non-trivial effects on the level of 
cell shape and forces. In order to demonstrate this important aspect, we have studied two 
important additional features of the cytoskeleton. First non-linear links were introduced 
via the worm-like-chain model (WLC). In the WLC case the passive square network first 
contracts and then expands again as tension is increased. The ACN requires much larger 
values of tension to contract compared to the linear case. Otherwise the arc morphology is 
essentially the same as in the linear case. 

As a biologically very relevant aspect of the cellular cytoskeleton, we also have studied the 
adaptation response of network links. For certain parameter values, a saturation response for 
both elasticity and tension leads to a strong difference between EA and T of boundary links, 
which are strongly increased, and EA and T of internal links, which are increased much less. 
Similar aspect have been addressed before in the framework of Finite Element Modeling 
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(FEM) |40|. In this case, the biochemical regulation has been modeled with more detail. 
Both the resulting cell shapes and the formation of stress fibers inside the cell demonstrate 
that the FEM model strongly resembles the HSN studied here. Therefore it would be 
interesting to combine the detailed biochemical model with the actively contracting cable 
network studied here. While it appears to be very challenging to develop a homogenization 
strategy for the ACN, it is interesting to consider if similar features as resulting from the 
ACN could be obtained in a FEM-framework. 

Active cable networks have also been shown to describe the circular arc morphology 
of tissues pinned at discrete sites [4^|. Because here arc radius also scales with spanning 
distance as for the arc radius of strongly adhering cells, the tension-elasticity model seems 
to capture all essential element of this situation. In the tissue case, the cable network 
represents the fibrous nature of the collagen matrix and the active contractility corresponds 
to cell contraction. Because in addition water can flow out of the contracting cell-matrix 
composite, volume is not conserved and compression is not propagated. Therefore the 
standard models of elasticity are strictly speaking not appropriate. Indeed they do not 
predict the circular arc morphology, but rather show flat contours corresponding to the 
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unconstrained reference shape of the elastic model 

In summary, HSN, PCN and ACN are simple model systems which however show sur- 
prisingly rich responses to internal contractility and therefore lead to interesting conclusions 
about the physical elements required to endow cells with a sense of geometry. ACN seem to 
be very appropriate to model strongly adherent tissue cells as they not only implement some 
of the most important fundamental features of the cytoskeleton (asymmetry under tension 
and compression, contraction by molecular motors), but also lead to functions which are 
very reminiscent of real cells (robustness under structural re-arrangements and adaptation 
to local adhesion constraints). 
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